load c:\manoj\projects\plasma\HUA_PIU1.mat
load c:\manoj\projects\ace\ESKLER fday_esk_ler_ts_filt FILT_LER FILT_ESK;
load c:\manoj\projects\plasma\Julia_PLP.mat Julia_W fday w;
load c:\manoj\projects\ace\OMNI_ELEC ace_all; % JUST ACE IEF
edit ace_julia_pcn.m
load c:\manoj\projects\ace\pcn; mean_delta = 6.9444e-004;
jfday = fday(1,:);
LL = magnetic_frac >= jfday(1)-floor(jfday(1)) & magnetic_frac <= jfday(end)-floor(jfday(end));
mjd_date = datenum(2000,1,1);
for i = 221:749,
    subplot(5,1,1);
    set(gca,'FontSize',16)
    jfday = fday(i,:);
    L = ace_all(:,1) > jfday(1) &  ace_all(:,1) < jfday(end);
    plot(ace_all(L,1),ace_all(L,2),'LineWidth',2);
    title('ACE Ey');
    ylabel('mV/Km');
    set(gca,'XTick',[])
    subplot(5,1,2);
    set(gca,'FontSize',16)
    L=PN_fday > jfday(1) & PN_fday < jfday(end);
    plot(PN_fday(L),PN(L),'LineWidth',2);
    title('PCn index');
    set(gca,'XTick',[])
    ylabel('nT');
    subplot(5,1,3);
    set(gca,'FontSize',16)
    L = fday_esk_ler_ts_filt > jfday(1) &  fday_esk_ler_ts_filt < jfday(end);
    plot(fday_esk_ler_ts_filt(L),FILT_LER(L),'LineWidth',2);
    title('LER magnetic data');
    ylabel('nT');
    set(gca,'XTick',[])
    subplot(5,1,4);
    set(gca,'FontSize',16)
    plot(jfday,w(i,:),'LineWidth',2);
    title('JULIA');
    set(gca,'XTick',[])
    ylabel('V_z (m/s)');
    subplot(515);
    set(gca,'FontSize',16)
    [a,b] = intersect(HUA_PIU_fday,floor(jfday(1)));
    if ~isempty(b),
    plot(HUA_PIU_fday(b)+magnetic_frac(LL)++mjd_date,HUAPIU(b,LL),'LineWidth',2);
    title('HUA_{\DeltaH}-PIU_{\DeltaH}');
    ylabel('nT');
    xlabel(datestr(floor(mjd_date+jfday(1))));
    end;
    
    datetick;
    pause;
end;